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No spin-singlet bb quarkonium state has yet been observed. In this paper 
we discuss the production of the singlet P-wave bb and cc 1 P\ states h b and 
h c . We consider two possibilities. In the first the 1 P± states are produced 
via the electromagnetic cascades T(3S) — > T]b(2S) +7 — > /ift + 77 — > 776 + 777 



. and ip' ~ ^ ^c+77 — ► ^c+777- A more promising process consists of 

single pion transition to the l Pi state followed by the radiative transition to 
the I 1 S state: T(3S) -> /i fo +7r° -> ?7 6 +7r +7 and?/'' -> /i c +vr° -> ?7 c +7r +7. 
For a million T(3S') or ^"s produced we expect these processes to produce 
several hundred events. 

a: 

PACS Categories: 14.40. Gx, 13.20.Gd, 13.40.Hq, 12.39.Ki 

<d : 

The study of bound states of heavy quarks has provided important tests of quantum 
chromodynamics (QCD) PJ. The heavy quarkonium cc and bb resonances have a rich 
spectroscopy with numerous narrow S, P, and D-wave levels below the production 
threshold of open charm and beauty mesons. The spin-triplet S'-wave states, ip(nS) and 
T(nS) with J PC = 1 , are readily produced by virtual photons in e + e _ or hadronic 
interactions, and then undergo electric dipole (El) transition to the spin-triplet P-wave 
levels. Previous studies have discussed the production of the spin-triplet D-wave bb 
states p], H and there has been some discussion of how one might produce the l*Pi cc 
state |5], |5], [7|, Up to now, the only observed heavy quarkonium spin-singlet 

state has been the ?7 c (l 1 S'o), but the Belle Collaboration |L0| has just announced the 
discovery of the ^(2 1 5'o) in B decays at a mass of (3654 ± 6 ± 8) MeV/c 2 . There 
have also been a few measurements suggesting the l 1 Pi(cc) state in pp annihilation 
experiments fll , |i~2l , p~3|| but these results have yet to be confirmed. No bb spin-singlet 



states have yet been seen. 

The mass predictions for the singlet states are an important test of QCD moti- 
vated potential models [0, |I§ [16], [T7|, 0, [[J, |20|, |2Tj, ^, g3| and the applicability 
of perturbative quantum chromodynamics to the heavy quarkonia cc and bb systems 
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!5|, E7|, as well as the more recent NRQCD P8|| approach. For QCD-motivated 
potential models the triplet-singlet splittings test the Lorentz nature of the confining 
potential with different combinations of Lorentz scalar, vector, etc., giving rise to differ- 
ent orderings of the triplet-singlet splittings in the heavy quarkonium P-wave mesons. 
Furthermore, the observation of cc and bb states and the measurement of their masses 



is an important validation of lattice QCD calculations [29 , [30| , |3T| , |32| , |33| , 34] , which 
will lead to greater confidence in their application in extracting electroweak quantities 
from hadronic processes. Under the assumption of a Fermi-Breit Hamiltonian and only 
vector-like and scalar-like components in the central potential, Stubbe and Martin |35] 



predicted that the n x Pi mass lies no lower than the spin-averaged 3 Pj masses (weighted 
with the factors 2J + 1), denoted by n 3 P cog . Violation of these bounds would indicate 
a significant underestimate by of relativistic effects. 



In Table | we summarize some predictions for hyperfine mass splittings for P-wave 
cc and bb levels. The wide variation in the predicted splittings demonstrates the need 
for experimental tests of the various calculational approaches. 

There are two possibilities for producing spin-singlet states. In the first, the system 
undergoes a magnetic dipole (Ml) transition from a spin-triplet state to a spin-singlet 
state. The predictions for Ml transitions from the T(n 3 5'i) levels to the rjb(n a So) 
states, for both favored Ml transitions and hindered Ml transitions with changes of 
the principal quantum number, have been reviewed in Ref. [Q. The second route 



begins with a hadronic transition, from a n 3 Si state to a 1 P\ state, emitting one or 
more pions, followed by the electromagnetic decay of the 1 P\ state. 

In this paper we examine the production of the spin-singlet P-wave cc and bb states. 
We examined the decay chains that start with the Ml transition from the ip' to the 
r]' c in the cc system and from the T(3S) to either the r] b (3S) or r] b (2S) state in the 
bb system. In both cases the Ml transition is followed by an El transition to the 
spin-singlet 2P or IP state. This is in turn followed by a second El transition to a 
t^Sq state. In addition, the 2 1 Pi bb state can undergo an El transition to the l 1 !^ 
state. However, with the current CLEO data set, the only decay chain which has any 
hope of being seen in the bb system is T(35) — > rjb(2S)j — ► hf,(lP)yy. We therefore 
only present results relevant to this set of decays. 

The decay chains originating with the hadronic transitions are more promising. We 
therefore include estimates of branching ratios for chains originating with the direct 



hadronic transition T(3S f ) — > h^Pi) + vr° discussed by Voloshin |J7| followed by the 
radiative decay hb( l Pi) — > r] b (lS) +7 and the analogous transitions in the charmonium 
system ip'{2S) — > /i c ( 1 Pi) + 7r° — > ■q c {lS) n f'K Q . Kuang and Yan |38j have also considered 



the related spin-flip transition T(3S f ) — > h^Pi) +7T7t which may provide an additional 
path to the h b . 

Searches for the 1 P\ states have taken on renewed interest because of the current 
data-taking runs of the CLEO Collaboration at the Cornell Electron Storage Ring 
(CESR), which are expected to significantly increase their sample of data at the T(35*) 
resonance, and the proposed CLEO-c project which will study physics in the charmo- 
nium system. 

We begin with the bb mesons and decay chains involving only radiative transitions. 
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Table I: Predictions for hyperfine splittings M(n 3 P cog ) — M(n l Pi) for cc and bb levels. 



Reference 


Approach n = 1 cc 


n = 1 bb 


n = 2 bb 




(MeV) 


(MeV) 


(MeV) 



GI85 g 



MR83 P] 
LPR92 P| 



OS82 
MB83 fll8 



GRR86 
1087 m 
GOS84 ?7 S = 1 
GOS84 rj s = H 
PJF92 IE3I 



21 



HOOS92 [|D 
PTN86 
PT88 



2| 



a. 
b 
c 
d 
c 

f 

g 
h 
h 

i 
j 
j 
J 




4 
10 

-5 
-2 
24.1 ±2.5 
6 

17 

-20.3 ±3.7 
-0.7 ±0.2 

-3.6 

-1.4 



2 


2 
3 

-2 
-1 
3.73 ±0.1 

3 

8 

-2.5 ± 1.6 
-0.18 ±0.03 

-0.4 

-0.5 



2 
1 
1 

3 

-2 
-1 

3.51 ±0.02 
2 
6 

-3.7 ±0.8 
-0.15 ±0.03 
-0.3 
-0.4 



SESAM98 31 


k 




~ -1 




CP-PACS00 [33 


1 


1.7-4.0 


1.6-5.0 





a Potential model with smeared short range hyperfine interaction. 
(The splittings are based on masses rounded to 1 MeV, not the results 
rounded to 10 MeV as given in Ref. [QO 

6 Potential model with long range longitudinal color electric field. 
c Potential model with PQCD corrections to short distance piece. 
d Potential model with smeared hyperfine interaction. 

e Potential model with smeared hyperfine interaction and relativistic corrections 
* Potential model includes 1-loop QCD corrections. 
9 Potential model with short distance from 2-loop PQCD calculation. 
Results shown for A Ms = 200 MeV. 

h Potential model with confining potential with both Lorentz scalar and vector. 
r) s gives the fraction of the confining potential that is pure Lorentz scalar 
versus Lorentz vector. 

1 Potential model. Solution is for Richardson potential and m c = 1.49 ±0.1 GeV. 



Other solutions given in Ref. [^2J are consistent with this result within errors. 
* PQCD 

k Unquenched nonrelativistic lattice QCD. 

1 Lattice QCD; the result is dependent on the value used for (3 and uiq. 
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Table II: Radiative electric dipole transitions involving hb(l l P\) and h' b (2 1 Pi) bb states. 
The details of the calculation are given in the text. 



Transition 


M< 


M f 


u 


(r) 


r 




(MeV) 


(MeV) 


(MeV) 


(GeV- 1 ) 


(keV) 


3% -> 2 l P l 


10337 


10258 


78.7 


-2.46 


3.2 


3 1 S -> l l Pi 


10337 


9898 


430 


0.126 


1.4 


2 1 P 1 -> l 1 ^ 


10258 


10148 


109 


-1.69 


2.7 


2 x Px -> 2% 


10258 


9996 


259 


1.57 


15.4 


2 1 P 1 - 1% 


10258 


9397 


825 


0.222 


10.0 


2% - l 1 ^ 


9996 


9898 


97.5 


-1.53 


2.3 


1 1 P 1 - 1% 


9898 


9397 


488 


0.940 


37 



To estimate the number of events expected from these decay chains we need to estimate 
the radiative partial decay widths between states and the hadronic partial widths of 
the appropriate 1 S and 1 P 1 states. 

The Ml transitions from the T(3S) to the rjb{3S) and rjb(2S) were studied in detail 



in Ref. [06 |, which we will use in what follows. The El transitions are straightforward 



to work out and in the nonrelativistic limit are given by 

rfSo - 'p, + 7 ) = ~« 4 ^ K^iMW (1) 

F( 1 P 1 -> % + 7 ) = ~« e% to 3 | (^rf^) | 2 (2) 

TCP, -+ 1 D 2 + 1 ) = % -a e% uj 3 \{ l D 2 \r\ l Pi)? (3) 

where a = 1/137.036 is the fine-structure constant, eg is the quark charge in units 
of |e| (—1/3 for Q = b), and u is the photon's energy. The photon energies, overlap 
integrals, and partial widths for the El transitions between l P\ and 1 5'o levels are given 
in Table [TT] and summarized in Fig. 1, along with the relevant Ml transitions. The n^S^ 
masses were obtained by subtracting the predictions of Ref. |14 for the n 3 S\ — ^Sq 
splittings from the measured n 3 Si masses, while the n 1 P 1 masses were obtained by 
subtracting the predictions for the n 3 P cog — n l Pi splittings of Ref. ]L4] from measured 



n 3 P cog values. The overlap integrals, (r) = ( 1 L^,|r| 1 L L ), were evaluated using the 
wavefunctions of Ref. |T3| . We found that the relativistic effects considered in Ref. |14 
reduce the partial widths by a few percent at most. Somewhat larger matrix elements 
were obtained in an inverse-scattering approach 0, except for the highly suppressed 
3S — > IP transition, whose matrix element is very sensitive to details of wave functions 

To estimate the number of events in a particular decay chain requires branching 
fractions which depend on knowing all important partial decay widths. Inclusive strong 
decays to gluon and quark final states generally make large contributions to the total 
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Figure 1: Radiative transitions in the bb system. The dashed lines represent Ml 
transitions, the solid lines El transitions and the dotted lines single 7r° emission. The 
transitions are labelled with their partial widths given in keV. 
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width and have been studied extensively HO, II , f|5 1 . The relevant theoretical 

expressions, including leading-order QCD corrections [41], are summarized in Ref. [44]: 



r( 1 s -^) = ^/W0)| 2 (4) 

with a multiplicative correction factor of (1 + 4.4^) for bb and (1 + 4.8—) for cc, 

rCPi - 999) = ^\R' P (0)\ 2 Hm Q (r}) (5) 

and 

T( 1 P 1 - ^ + 7 ) = - eJ-r( 1 P 1 - ^) (6) 

5 y a s 

where we also include the decay 1 Pi — ► (7(7 + 7. 

Considerable uncertainties arise in these expressions from the model-dependence of 



the wavefunctions and possible relativistic contributions |14| . In addition, the loga- 
rithm in the decay T( 1 Pi — > ggg) is a measure of the virtuality of the quark emitting 
the gluon. Different choices have been proposed for its argument, introducing further 
uncertainty. Rather than evaluating these expressions in a specific potential model, 
we can obtain less model-dependent estimates of strong decays by relating ratios of 
theoretical predictions, in which much of the theoretical uncertainties factor out, to 
experimentally measured widths. Although we expect the wavefunction at the origin 
to be slightly larger for the singlet state than the triplet state, we expect this difference 
to be much smaller than the uncertainties mentioned above. 

To make our estimates, we will need in addition to Eqs. (3)-(5), the following 

n 3 s 1 ^ 999 ) = m ^ )a h m\ 2 (7) 



expressions 



with a multiplicative correction factor of (1 — 4.9—) for bb and (1 — 3.7—) for cc, 



n i s 1 7 + gg) = y n ^ s \m\ 



, _ 32(^-9)4^, 

9m3 



with a multiplicative correction factor of (1 — 7.4 2a ) for bb and (1 — 6.7 s *) for cc, and 

F( 3 Pi ^qq + g) = ^|P'p(0)| 2 ln(m Q (r)) , (9) 

Q 

where the QCD correction factor for the last expression is not known. Taking account 
of decays of the 2 3 Si and 3 3 5*i states to 7i7iT(nS), lepton pairs, and Xb(nP)j, as quoted 
in Ref. |j46| , we find total branching ratios to non-glue final states, and assume glue 
to constitute the remainder. Using branching ratios and total widths quoted in Ref. 
|f46| , we then arrive at the estimates summarized in Table |TT] for r[T(nS) — ► glue] = 
r[T (nS) — > hadrons] + T[T(nS) — > 7 + hadrons]. 



6 



Table III: Ingredients in estimates of T[T(nS) —>■ glue]. 



T(nS) 


£>(non — glue) 


B(glue) 


r(tot) 


r(glue) 


state 


(%) 


(%) 


(keV) 


(keV) 


T(2S) 


49.1 ± 1.5 


50.9 ± 1.5 


44 ±7 


22.4 ±3.6 


T(3S) 


44.9 ± 1.4 


55.1 ± 1.4 


26.3 ±3.5 


14.5 ±2.0 



Using Eqs. (7) and (8), a s (T 2S ) = 0.181, and a s (T 3S ) = 0.180 we find r[T(25) -> 
hadrons)] = 21.8 ± 3.5 keV and T[T(35) -> hadrons)] = 14.1 ± 2.0 keV. The ratio of 
the widths from Eqs. (4) and (7): 

27-7T 1 (4 + 4 4^ 

results in T[r] b (2S) -> hadrons)] = 4.1 ± 0.7 MeV and T[r) b (3S) -> hadrons)] = 2.7 ± 
0.4 MeV. 

We follow the same procedure to estimate the hadronic width for the l P\ states, 
although in this case we need to make use of a theoretical estimate for the partial width 
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Pi — > Si 7. Here we have |fP 



T( 1 P l -> hadrons) = — x T( 3 Pi -> hadrons) (11) 

where n./ is the number of light quark flavours in the final state which we will take 
to be 3, ignoring the kinematically suppressed charm-anticharm channel. This results 
in a conservative upper limit for T( 1 Pi — > hadrons) and hence a lower limit for the 
branching ratio of this state to 7 ± rj b . As mentioned, the QCD corrections to these 
widths are not known. The large uncertainties arising from the wavefunction and 
logarithms in Eqs. (5) and (9) cancel out. 

The only branching ratios quoted in Ref. [IB] for the n 3 Pi states are for decays 



to n /3 Si7. Using quark model predictions for the radiative transitions and assuming 
that hadronic decays dominate the remainder of the total widths, we can estimate the 
hadronic partial widths of these states. The results are summarized in Table |TV| . 

The branching ratios obtained by combining the partial widths given in Tables [II 
and [TV] are summarized in Table |V[ 

To study the singlet P-wave bb states we considered the two-photon inclusive tran- 
sitions 3 3 Si ^ 3% ^ 2 X P X or ^ 1 1 P 1 and 3 3 Si ^ 2 1 S ^ l 1 ^. In all cases the x Pi 
states can undergo further PI radiative transitions to 1 S'o states. It may be that this 
last photon provides a useful tag to distinguish the cascade of interest from other pos- 
sible decays involving triplet P and P-wave bb states. We use the branching ratios pre- 
dicted in Ref. |6| for the initial Ml transitions, B(T(3S) -> rj b (3S) ±7) = 0.10 x 10~ 4 
and B(T(3S) — > r]b(2S) + 7) = 4.7 x 10~ 4 , which correspond to the GI mass-splittings 
and wavefunctions and where the latter result takes into account relativistic cor- 
rections. Combined with the branching ratios for the subsequent PI transitions given 
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Table IV: Partial widths of 3 Pi and l P\ states. The details of the calculation are given 
in the text. 



3 P 


E„£( 3 Pi^n 3 S l7 ) a 


r( 3 Pl _> s Sil) 


r(— > hadrons) 


state 


(%) 


(keV) 


3 P 1 (keV) 1 P 1 (keV) 


l 3 Pi(66) 


35 ±8 


32.8 b 


60.9 50.8 






28.9 C 


53.7 44.7 


2 3 P 1 (bb) 


29.5 ±4.2 


25.2 6 


60.2 50.2 






16.8 C 


40.1 33.4 



a Particle Data Group [4§ 
b Kwong and Rosner 



Godfrey and Isgur [14] 



Table V: Partial widths and branching ratios for spin-singlet bb states. The details of 
the calculation are given in the text. 



Initial 


Final 


Width 


B 


state 


state 


(keV) 


(%) 


3% 


2 x Pi7 


3.2 


0.12 




l X Pi7 


1.4 


0.05 




99 


2700 


99.8 


2 1 P 1 


2% 7 


15.4 


19.3 




1^07 


10.0 


12.5 






2.7 


3.4 




999 


50.2 a 


62.8 




199 


1.6 


2.0 


2 1 5' 


l 1 P ll 


2.3 


0.057 




99 


4100 


99.9 


l'Pi 


1%7 


37.0 


41.4 




999 


50.8 a 


56.8 




199 


1.6 


1.8 



a Based on the partial width for 3 P 1 -> 3 S^ of Ref. @ in Table [TV]. 
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in Table V the only decay chain that might yield enough events to be observed is 
3 3 Si 2 1 Sq l x Pi which yields roughly 0.3 events per million T(3S) states pro- 
duced. 

A more promising approach is the decay chain T(3S) — > 1 P 1 7r° followed by the 
El radiative transition 1 P 1 — > 1 S'o7. Voloshin estimates B(T(3S) — > l x Pi + 7r°) = 
0.10 x 10~ 2 |7|. Thus, fi[T(3S) -> l l Pi + tt° -> 1% 7 ] ~ 4 x 10~ 4 , which would yield 



~ 400 events per million T(35) produced. This signature should be easily seen by the 
CLEO detector, which has excellent photon detection capabilities. Since the recoil of 
the l x Pi state is relatively small, the 488 MeV photon from the l l P\ — > l 1 ^ decay 
(suitably Doppler-shifted by up to ± 20 MeV) should provide a useful tag. Kuang and 
Yan predict |38| the partial width for the hadronic transition T(3S) — > hb(l l Pi) + tttt 



to be 0.1-0.2 keV, giving a branching ratio of ~ (3.8-7.6) x 10~ 3 . This is substantially 
higher than the value for B(T(3S) — > l x Pi + ir°) quoted above so it could provide an 
alternative path to the hj,- However, Voloshin [R7| does not obtain such a favorable 
branching ratio for this process, finding instead < 10 -4 . 

We now turn to the charmonium system. The search for the h c was discussed 
recently by Kuang || so we will be brief in our analysis, emphasizing aspects that are 
different from Kuang ||. As in the case of bb there are two routes to the h c . The first 
is the decay chain if>' — > 77^,7 — > and the second is through the hadronic transition 
V>' -> h c n°. 

For the first case we need the various radiative widths. The expression for the El 
width is given by Eq. (1), while the rates for magnetic dipole transitions are given in 
the nonrelativistic approximation by 

A 2 

r( 3 ^ % + 7 ) = ^a; 3 |(/|io(W2)N)| 2 (12) 

A 2 

ifSo - 3 ^ + 7) = ^uj 3 \(f\Ukr/2M 2 (13) 

m Q 

where we take m c = 1.628 GeV. The results, using the wavefunctions and l x Pi mass 
of Ref. are summarized in Table |VT[ To calculate T(ip' —>■ rj^rf), we took M(j]' c ) = 



3654 MeV, the central value quoted in Ref. ]Ll| ■ Note that the widths for the hindered 
Ml transitions are very sensitive to the wave functions. The hadronic widths for the 
7]' c and h c given in Table [VT] were obtained using the same procedure used for the bb 
hadronic widths: We relate theoretical expressions for ratios of the widths to a known 
measured width and take a s (ip') = 0.236. (In contrast to the bb system, the total width 
of the l 3 Pi cc meson is known J46|: r tot (x cl ) = 0.88 ±0.14 MeV.) The predicted result 
for h c — > hadrons is consistent with the NRQCD result obtained by Bodwin, Braaten 



and Lepage ||7|]. Combining these results we find that B(ip' — > rj' c j) x B{rj c — > h c ^) ~ 
10 -6 , which would yield a modest number of h c mesons at best. 

As in the case of the h^, a more promising avenue is the single pion transition 
1/}' — > h c 7r° followed by the radiative transition h c — > r) c "f, where the photon is expected 
to have an energy very close to 496 MeV and can be used to tag the event. Using the 
branching ratio of B(if)' — > h c ir°) = 0.1% predicted by Voloshin ]37j (see also Ref. [ff|,|38||) 
and the branching ratio given for h c — > rjc'y in Table [VT] , we obtain B(ip' — > h c -K°) x 
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Table VI: Partial widths and branching ratios for spin-singlet cc states. In column 5, 
O represents the operator relevant to the particular electromagnetic transition; O = r 
(GeV -1 ) for El transitions and O = jo(kr/2) for Ml transitions. The details of the 
calculation are given in the text. 



Initial 


Final 


Mi 


M f 


to 


(f\o\i) 


Width 


B 


state 


state 


(MeV) 


(MeV) 


(MeV) 




(keV) 


(%) 


2 3 S 1 


2% 7 


3686 


3654 


31.8 


0.982 


0.051 


0.018 




1%7 


3686 


2980 


638 


0.151 


9.7 


3.5 


2% 


HP l7 


3654 


3517 


134.4 


-2.21 


51.3 


0.69 




l 3 S l7 


3654 


3097 


515 


-0.0973 


6.3 


0.084 




99 










7400 


99.2 




1%7 


3517 


2980 


496 


1.42 


354 


37.7 




999 










533 


56.8 




199 










52 


5.5 



B(h c — > r] c j) = 3.8 x 10~ 4 , which is substantially larger than the decay chain proceeding 
only via radiative transitions. In his recent paper Kuang |J finds h c production to be 
sensitive to 3 Si — 3 Di mixing, so that a measurement of B(ip' — > h c ir ) would be a 
useful test of detailed mixing schemes between the ip' and the ^(3770) = ip", some of 
which are discussed in Ref. [48. 49 . 



Another promising approach for the detection of the h c has recently been proposed 
by Suzuki ||. He suggests looking for the h c by measuring the final state •yr] c of 
the cascade B — > h c K/K* — > / -fi] c K/K*. This channel is especially timely given the 
announcement by the Belle Collaboration of the discovery of the rj' c (2 1 So) in B decays 



1(J and, previously, the observation of the related decay, B — > xoK [0 • 

In the case of the S-wave ( 1 5'o) states, one should also bear in mind that 77 collisions 
have been used to observed the r\ c in several experiments (see |46|]). One candidate for 
77 - > Vb with mass 9.30 ± 0.02 ± 0.02 GeV/c 2 (consistent, however, with background) 
has been reported by the ALEPH Collaboration |Hfl . 

To conclude, we have explored different means of looking for the l P\ states in heavy 
quarkonium. In both the bb and cc systems the l x Pi state can be reached via the chain 
3 Si — » 1 S + 7 — > l 1 Pi + 77. However, in both systems one only expects of the order 
of a few events per million T(35) or ip' produced. In both systems, a more promising 
avenue is the transition 3 S\ — > 1 P\ + 7r° followed by the El radiative transition to the 
l 1 ^ state which would yield several hundred events per million T(3S") or ^"s produced. 
The alternative suggestion |38j of searching for the transitions 3 Si — > 1 P\ + tttt also is 
worth pursuing. 

One of us (J.L.R.) thanks Steve Olsen and San Fu Tuan for informative discus- 
sions. The authors thank Christine Davies for helpful communications. This work 
was supported in part by the United States Department of Energy through Grant No. 
DE FG02 90ER40560 and the Natural Sciences and Engineering Research Council of 
Canada. 
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